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Abstract 

The connection between the enhancement factor + of the pho- 
tonuclear El sum rule and the orbital angular momentum g-factor (gi) 
of a bound nucleon is investigated in the framework of the Landau- 
Migdal theory for isospin asymmetric nuclear matter. Special em- 
phasis is put on the role of gauge invariance to establish the n — gn 
relation. By identifying the physical processes which are taken into 
account in k and g^, the validity and limitations of this relation is dis- 
cussed. The connections to the collective excitations and to nuclear 
Compton scattering are also shown. 
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1 Introduction 



The enhancement factor (1 + k,) of the photonuclear El sum rule and the 
orbital angular momentum g-factor (g^) of a bound nucleon have attracted 
the attention of nuclear theorists as well as experimentalists for a long time, 
since these quantities reflect the presence of exchange forces and mesonic 
degrees of freedom in nuclei [1, 2, 3, 4]. More than 30 years ago, Fujita 
and Hirata [5] used the isospin symmetric Fermi gas as a model for an N=Z 
nucleus to derive the simple relation 1 + k = 2 g^^iy between k and the 
isovcctor (IV) part of in first order perturbation theory. Later it has been 
shown [6] that, because of the presence of correlations between the nucleons, 
only a part of the total k is related to ge^iy- It has been argued [7] that this 
part of K is related to the sum of the El strength in the region of the isovector 
giant dipole resonance (GDR). In more recent years [8], this modified k — ge 
relation has been used to analyse the results of photo-neutron experiments 
[9] and photon scattering experiments [10, 11], in particular for nuclei in the 
lead region. In these analyses, the corrections to ge^iy were associated with 
the meson exchange currents. 

On the other hand, as early as 1965, Migdal and collaborators [12] used 
an approach based on a a gas of quasiparticles to relate k to the parameters 
characterizing the interaction between the quasiparticles (the Landau-Migdal 
parameters) . Combining this relation with the more general one between g£^iy 
and the Landau-Migdal parameters [13], their approach suggested that the 
relation 1 + k = 2 guy holds more generally without recourse to perturba- 
tion theory. The fact that their result involves the total k instead of just a 
part of it reflects the quasiparticle gas approximation. Concerning the or- 
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bital g-factor, however, the Landau-Migdal approach suggests that it is the 
total g^^iv, and not only the part associated with the meson exchange cur- 
rents, which enters in the relation to the El enhancement factor. By using a 
somewhat different approach, this result has also been obtained by Ichimura 
[14]. 

The main advantage of the Landau-Migdal theory [13], which is based on 
the Fermi liquid approach due to Landau [15, 16], is that symmetries, like 
gauge invariance and Galilei invariancc, arc incorporated rigorously, and that 
the description of the collective excitations of the system is physically very 
appealing. The Fermi liquid approach to discuss sum rules in nuclear matter 
has therefore turned out to be very fruitful, and has been used in several 
papers on giant resonances [17, 18]. However, to the best of our knowledge, 
a general discussion of the physical processes which contribute to the K-g£ 
relation, as well as a discussion on the physical origin of this relation, is still 
outstanding. 

The strong interest in nuclear giant resonances and sum rules is continuing 
nowadays [19, 20], and new phenomena like the double giant resonances [21] 
or dipole resonances in neutron rich nuclei [22] have attracted attention. 
The Landau-Migdal theory, which is a strong candidate to describe the giant 
resonances [23] , is now extended in various directions [24] so as to give a more 
general description which is valid up to higher excitation energies. In the light 
of these recent developments and the analyses of photonuclear experiments 
mentioned above, more detailed understanding on the K-g^ relation would be 
desirable. 

In this paper we will present a general discussion on the K-g^ relation in 
isospin asymmetric nuclear matter, using the language of the Landau-Migdal 
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theory. The aims of our work are as follows: First, we will extend the relations 
obtained previously for the orbital g-factor and the El enhancement factor to 
the case oi N ^ putting special emphasis on the role of gauge invariance. 
Second, we will identify the physical processes which are taken into account 
in the K-g^ relation, both in terms of Feynman diagrams as well as time- 
ordered Bethe-Goldstone diagrams. In this connection we will also discuss 
the relation to the collective excitations. Third, we will show the connection 
to the photon scattering amplitude, which establishes the physical origin of 
the K-g^ relation. 

The rest of the paper is organized as follows: In sect. 2 we briefly re- 
view some relations of the Landau-Migdal theory which will be used for the 
discussions on gi^ and k. In sect. 3 we derive the expressions for the proton 
and neutron orbital g-factors in isospin asymmetric nuclear matter, which 
are in principle exact and hold also in relativistic field theory. In sect. 4 
we discuss the El sum rule and its relation to the orbital g-factors, limiting 
ourselves to the nonrelativistic case because of the problems associated with 
the center of mass motion. In sect. 5 we discuss the physical processes which 
are taken into account in the k — gi relation and establish the connection to 
the collective excitations. In sect. 6 we discuss the relation to the nuclear 
Compton scattering amplitude, which provides the physical origin of the K-g^ 
relation. A summary and conclusions are given in sect. 7. 
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2 Vertices and correlation functions in the 
Landau-Migdal theory 



In this section we briefly review some relations of the Landau-Migdal theory 
which will be used in later sections. 




Figure 1: Graphical representation of relation (2.1) for the vertex. The 
shaded square stands for the full vertex (F), and the black circle for the 
vertex which is two-body irreducible in the p-h channel (Fq). T denotes the 
T-matrix. 

The integral equation for the electromagnetic vertex F is graphically 
shown in Fig. 1 and is written symbolically ^ as 

V = To-tTSSVo, (2.1) 

where Fq is the part which is two-body irreducible in the particle-hole (p-h) 
channel, S is the full nucleon propagator, and T is the two-body T-matrix 
which satisfies the Bethe-Salpeter equation 

T^K-iTSSK (2.2) 

with the two-body irreducible kernel K. The important step in the renor- 
malization procedure of the Landau-Migdal theory is to split the product S S 

^In the symbolic notations we will omit the Lorentz indices for the electromagnetic 
vertices, currents and correlation functions. 
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into the two parts [13] 

SS^A + B, (2.3) 

where A denotes the product of the pole parts of the particle and hole prop- 
agators {A = ^^P"''^)^''"''^), 

including the prescription to evaluate the other 
energy dependent quantities in a frequency loop integral for on-shell p-h 
states. (The explicit form of A in the momentum representation will be 
specified below.) All other parts, like non-pole parts, the product of two par- 
ticle or two hole propagators, antinucleon propagators etc, are included in 
the quantity B. The basic idea here is that the part A represents the "active 
space" , while the effects of B are renormahzed into the effective vertex and 
effective interaction. 

The equation (2.2) for the T-matrix is then equivalent to 

T = T^^) - iT AT^""^ (2.4) 
tH = K-iT^^^BK, (2.5) 

and equation (2.1) for the vertex is equivalent to 

r = r('^) -ir('")^r (2.6) 

r('^) = To -iT^'^^STo. (2.7) 

As is clear from these equations, the quantities T^^^ and T^^^ do not involve 
the product S^°^'^S^°^'^ in the intermediate states, that is, they have no p-h 
cuts, and eqs. (2.4) and (2.6) can be considered as RPA-type equations with 
T^^^ and T^^^ playing the role of the effective interaction and the effective 
vertex, respectively, acting in the space of p-h states. 

Another quantity which is of particular interest for the collective excita- 
tions of the system is the electromagnetic current-current correlation func- 
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Figure 2: Graphical representation of eq.(2.8). For explanation of the sym- 
bols, see the caption to Fig.l. 

tion, or shortly the correlation function, which can be expressed in coordinate 
space as follows (see Fig. 2) ^: 

i (0|T (j^(x')j^(x)) |0) = I Tr [r^(|/, x\ y') S{y\ z') V^{z\ x, z) S{z, y) ] , 

(2.8) 

or symbolically in the form 

n^iV^SST (2.9) 

The trace (Tr) in (2.8) stands for an integral over the space-time positions 
z' , y', y as well as the trace over spin-isospin indices. 

If we use eqs. (2.3), (2.6) and (2.7) in eq.(2.9), we obtain the separation 
of the correlation function into a part 11^, which involves p-h cuts, and a 
part n^, which involves the rest, like 2p — 2h cuts etc.: 

U^ir'-'^^ AT + iToBT^'^^ = Ua + Ub. (2.10) 

The use of eq.(2.6) for the vertex F then generates the RPA series for the 
part n^. 

^We note that, although this diagram is drawn in the p-h channel, no particular time 
ordering has been chosen, e.g; the product SS in (2.8) can also involve two particle 
(forward propagating) Green functions. 
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For illustration, we represent Ha by Fig. 3, and show some selected 
examples ^ for diagrams which contribute to Ua in Fig. 4. 




Figure 3: Graphical representation of the part in eq.(2.10). The shaded 
square represents F, the black square F^*^), and p (h) denotes the pole part 
of the particle (hole) propagator. 

The common feature of these diagrams is that they involve at least one 
electromagnetic vertex where a p-h pair is created by the external field. This 
is in contrast to the diagrams contributing exclusively to n^, for which we 
show some examples in Fig. 5. 

Wc now specify the form of A in momentum space, which corresponds to 
the product of the pole parts of particle and hole propagators: 

sr'\k) = , sr'\k) = TT^^ , (2.11) 

where and Zk are the quasiparticle poles and residues, and is the Fermi 
distribution function. From the structure of the RPA-type equations (2.4) 
and (2.6) it is clear that each loop integral is associated with one factor A. 

'^Thc diagrams of Figs. 4 and 5 should be iindca'stood as examples for timc-ordcrcd dia- 
grams in the context of a mean field approximation, like the Hartree-Fock approximation, 
for the single particle Green functions. As we will discuss in detail in sect. 5, the 2p-2h 
diagrams shown here give contributions to both Ha and Hb- Examples for self energy 
corrections are not shown in this figure, but are also discussed in sect. 5. 
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Figure 4: Some selected time-ordered diagrams contributing to Ha- The first 
line shows the non-interacting part and a meson-exchange current contribu- 
tion in first order perturbation theory, the second line shows some graphs 
which arise from 2p-2h states in the vertex F^'^^ and the third line shows 
some examples reflecting the 2p-2h states in the interaction T^^^ . The vertex 
where the external field couples to a dashed line stands symbolically for the 
two-body exchange currents. Self energy corrections also contribute, but are 
not shown in this figure. 



We therefore consider the following expression to represent a generic loop 
integral which involves a p-h cut: 

X F(^ + |,A;-|), (2.12) 
where F does not contain A, and = (a;, q) is the momentum transferred 
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Figure 5: Some selected time-ordered diagrams which contribute exclusively 

to n^. 



by the external field. By closing the ko contour in the upper plane, we get 



I 



+ 



nfe+(l - Uk-) 



—F (efc+ k+ ; (-a; + ek+) fc-) ] +1f , 



(2.13) 



u - {ek+ - Sk-) + i-f' 
where 7 = sign a;, Ip denotes the contributions from the poles of F in the up- 
per plane, and all quantities with subindex k± refer to the three momentum 
k± — k± q/2. We expand both functions F in (2.13) around u; — €k+ — efc_, 
and take only the first term in this expansion. This, together with the pre- 
scription to omit the term Xp, defines the quantity A{k, q): 

d^k {nk+ — Uk-) Zk+Zk- „, J I \ /oi/tN 

/ TTTT^ 7 T— — F{ek+ k+ ; €k- fc_ 2.14 



J 



(2.15) 



From this identification we obtain 



A{k, q) ^ -2ni5 ( ko - ^^4^) "^-^ P{.) , (2.16) 

2 J uj- {ek+ - Sk-) + 1-1 
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where P{uj) denotes the prescription to evaluate F at u = ek+—ek-- Together 
with the delta function in (2.16) this means that A effectively replaces the 
function F{k + - , /c — -) by its "on-shell" (o.s.) value defined by 

fco + ^°=efe+, A;o-^°=efe_. (2.17) 

Therefore, the contribution of a Feynman diagram, which has p-h cuts, to 
is obtained by replacing the energy dependent vertices and T-matrices 
by their on-shell values. The rest is contained in Es. (Later we will specify 
what this means in terms of time-ordered Bethe-Goldstone diagrams, and 
make contact to the examples shown in Fig. 4.) 

A similar form of A can be derived also in finite systems. In this case 
one uses the "A-representation" [13] instead of the momentum representa- 
tion, which involves a set of single particle wave functions ^x{E) chosen so 
as to diagonalize the single particle Green function for fixed frequency E: 
Gyx{E) = Sx>xGx{E). The form of A then becomes 

Ayxiko, u;) = -2niS (ko - ^) inynx)Z,Zx p 
\ 2 J uj- {ey - ex) + «7 

where P{oj) implies to evaluate the function F, which appears in the same 
loop as A, 8ki uj — ey — ex- 

Finally in this section, let us use the form (2.16) to write down the expres- 
sion for and the equation for F in momentum space. For this purpose, 
we introduce the quasiparticle current j^'^^ and the quasiparticle interaction 
by 

j^''\k+,k^) = ^Zk+Zk- F(-)(e,+ fc+; e,_ fe_) (2.19) 
&\k+,k_;e+,e^) = ^JZk+Zk-Zt^Z,_ T^^\ek+k+, efe_fc_; e,+£+, e,_^_) . 

(2.20) 
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Defining also the total current j(A;+, k-) in terms of the total vertex F in a 
way similar to (2.19), the expression for Ua can be written as 

nA(^,q) = v f ^tr U-\k.,k^) — ■ Jik^.k.)] 

J (27r)'^ V (J - {€k+ - €k-) + ^7 / 

(2.21) 

where V is the volume of the system, and the equation for the current takes 
the form 



j{k+,k_)=j^-\k^,k^) 



(27r) 



rtr 



(2.22) 



Relations analogous to (2.21) and (2.22) can also be written down for finite 
systems in the A-representation. 

3 The orbital g-factor 

In this section we will extend the previously obtained relations [13, 26] for the 
orbital g-factor of a quasiparticle on the Fermi surface to the case of isospin 
asymmetric nuclear matter [N ^ Z) . Although our notation will refer to the 
nonrelativistic case, the result is valid also in relativistic theories. 

The magnetic moment associated with the quasiparticle current j^'^^ is 
obtained from m = — | Vg x (^j^'^^S^^\p + q — p')) _o" nuclear matter, 
the only way to get a term proportional to the orbital angular momentum 
is to let the derivative act on the momentum conserving delta function [26]. 
Since the matrix element of the orbital angular momentum operator between 
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plane wave states is ipx 'V p6^^\p — p') , we obtain for the angular momentum 
g-factor of a proton or neutron at the Fermi surface 

9e{(^) = —r^ja^ (« =P:n), (3.1) 

where Cp > is the proton electric charge, Pf{p) and ppiji) denote the 
Fermi momenta of protons and neutrons, and the current refers to the 
g — > limit of the quasiparticle current (2.19) at the Fermi surface^. The 
calculation of the angular momentum g-factor in nuclear matter therefore 
reduces to the calculation of the quasiparticle current at the Fermi surface 
at zero momentum transfer. 

Invariance with respect to local gauge transformations of the charged 
particle fields leads to the Ward-Takahashi identity for the electromagnetic 
vertex, which in turn gives expressions for the currents ja in the limit — > 
followed by |qf| 0, which is called the "static limit", following ref. [13]. (The 
use of the Ward-Takahashi identity in many-body systems has been explained 
in detail in refs. [26, 27] for relativistic theories, but for convenience we 
summarize some relations in Appendix A.) The results for the static limit of 
the proton and neutron currents are 

if) = e,^^(p), jf^ = 0. (3.2) 

We denote by vf{ol) the Fermi velocities of protons {a — p) and neutrons 
{a — n). 

Since we need the quasiparticle currents j'^^ instead of the static currents 

*Here and in the following sections, currents without explicit momentum variables refer 

to the limit of zero momentum transfer. Also, in order to avoid unnecessary indices, we 
will frequently use V =V^ and T = T^^ for any vector V or tensor T^^ which depend on 
the momentum of the quasiparticle, which is put along the 3-axis. 
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we have to use eq.(2.22) in the hmit^a; — > followed by \q\ — > 0. In this 
limit, we can use 

^-Sm-PF) (3.3) 



for either protons or neutrons. Then both vectors k and t in eq.(2.22) are on 

the Fermi surface, and the integral reduces to an integral over the angle be- 
tween the directions k and i (Landau angle). In this kinematics, the interac- 
tion i^'^) depends only on the Landau angle and becomes the "Landau-Migdal 
interaction" . It is convenient to introduce the dimensionless interaction 

T{k-l) = '^—&\k-l), (3.4) 

where by definition we identify pp and vp with the Fermi momentum and the 
Fermi velocity in symmetric nuclear matter with the same baryon density. 
It is then easy to obtain the following relations from eq.(2.22) and (3.2) (see 
Appendix A): 

2 



Jp "^p 



MP) + f (^1 F,(pp) 



(3.6) 



where Fi{pp) and Fi{pn) — Fi{np) are the coefficients of the i — 1 Legendre 
polynomials Pi{k ■ £) in an expansion of the spin-independent part of Tpp or 

~pn 



Tpn ^. Lorentz invariance, or Galilei invariance in the nonrelativistic case. 



^Since the quasiparticle current j^'^^ by definition has no p-h cuts, the limits w — > 
followed by |q| commute for j'^^^ On the other hand, the total current including p-h 
cuts involves the quantity A of eq.(2.16), which vanishes if \q\ first, but assumes a 
nonzero value according to eq. (3.3) if a; — > first. 

^The relations to the more familiar parameters -Fi and F{ are Fi{pp) = Fi + and 



14 



can then be used (see Appendix A) to rewrite these expressions in the form 



Jp 



Pf(p) _ P£M}!LFJpn)^ ( ^^^^^ 



Pf 



6pF 



Pf 



(3.7) 
(3.8) 



where /ip and fin are the chemical potentials (Fermi energies including the 
rest masses) of protons and neutrons, and 



1 - 



N - Z 
A 



PF{p)pF{n) 

The orbital g-factors then can be expressed as follows ^ : 

M 



(3.9) 



Ap) = 

in) = 



Mvf 



Fi pn -/3 

ZpF l^p A 



(3.10) 
(3.11) 



We note that the following relation holds between the proton and neutron 
orbital g-factors: 



Z flp N , . Z 



(3.12) 



In a nonrelatistic theory we have f^p — fJ^n — M, and eq. (3.12) is the 

extension to N Z oi the well known fact [26, 27] that the isoscalar orbital 

g-factor is renormalized exclusively by relativistic effects. 

Besides the dependence on the neutron excess shown explicitly in (3.10) 

and (3.11), the Landau-Migdal parameter Fi{pn) might also depend on {N — 

^Eqs. (3.10), (3.11) are exact in nuclear matter, and include mesonic, relativistic, and 
configuration mixing effects. Because the vertex F^'^' is an effective quantity in the p-h 
space, conventional RPA-typc contributions, like the first order configuration mixing, are 
not included in (3.10) and (3.11). 
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Z^l A for fixed A. If the range of the interaction is small compared to 
this dependence can be expected to be weak. For the case of the long range 
1-pion exchange interaction, we note that in Fi{pn) the average over the 
cosine of the Landau angle (cosO^, = x) is performed with respect to the 
momentum transfer in the exchange channel, which is — (fc — = 



Prip) + Ppiny - 2pF{p)pF{n)x = 2pj,{l - x) + O 



( N-Z \ 

\ ) 



Therefore, 



up to O 4— j, the dependence on the neutron excess is as shown by the 
explicit factors N/A and Z/A in the relations (3.10) and (3.11). 

4 The El sum rule and the k — ge relation 

While the treatment of the orbital g-factor in the previous section was com- 
pletely general, the discussion of the El sum rule is plagued by problems 
related to the spurious center of mass motion [25]. Therefore, from now we 
have to restrict ourselves to a nonrelativistic framework, and consider the 
case of a Hamiltonian 



H 



2M 



\ f dV V^(r) A^(r) (4.1) 
A J 



+ IJ dVi J dV2 ^n^i)^n^2) ^(ri, ra) ^(r2)^(ri) (4.2) 
= Ho + V (4.3) 

with the two-body potential 

V{ri, rs) = Voiri, ra) + Ti • T2 Vi{ri, ra). (4.4) 

4.1 Effective charges and gauge invariance 

In this subsection, we digress to a discussion of local gauge invariance of a 
Hamiltonian containing two-body potentials. Our aims are, first, to investi- 
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gate whether the introduction of the famihar El effective charges, which ehm- 
inate the spurious center of mass motion, is consistent with gauge invariance, 
and, second, derive some relations for the 2-body part of the electromagnetic 
interaction which will be used later. 

Let us consider local gauge transformations of the form 



■0(r') — > exp iOxif) '0(^) 



(4.5) 



where x(r) is the gauge function, and 



1 + r. 



■Qp + 



l-r. 



Qr 



In the first place, of course, we take the physical electric charges 



(4.6) 



-e > 



and e„ = for Qp and Q„, but wc wish to keep our expressions more general in 
order to incorporate the effective charges in the later discussions. This gauge 
transformation effectively changes the first quantized operators as follows: 



exp 
exp 



-i {Q^hi + Q^'hj)] v{ij) [i (Q«xi + Q^^hj) 



(4.7) 



(4.^ 



where Xi = x(^i)) Xj = x{fj)j the other notations are self evident. 
Although it is not difficult to work out these expressions for a general gauge 
function, we consider here the long wave length limit (LWL) limit in order 

to derive the necessary relations most directly. In this case it is sufficient^ to 
consider x(r) = cq + r ■ c, where the constants. Using the fact that 

the 2-body potential commutes with the total 2-body isospin operator, the 



^This is easily seen [3] by noting that the vector potential in the LWL is a constant, 
and a gauge transformation A ^ A + Vx should again give a field of the same kind. 
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quantity [Q^^^Xi + Q^'^^Xj) the gauge transformation (4.8) can be replaced 
by c- (-Dq + Dq^^, where Dq — Q^'^^ri is the dipole operator for particle i, 
referring to arbitrary charges Qp and Qn- Since also (4.7) can be expressed 
by the dipole operator for particle i, the gauge transformation of the total 
Hamiltonian in first quantization 

H ^ Ho + V = Y,Ho{i) +Y.y{h j) (4.9) 

i i<j 

can be written in terms of the dipole operator Dq — ^ £>q as follows: 

i 

H exp [-ic ■ Dq] H exp [ic ■ Dq] . (4.10) 



This expression can be expanded in a series of multi-commutators: 



H^H-ic- [Dq, H] - \}(^ 



+ ... 



(4.11) 



where we have shown the terms up to C^(e^) explicitly. The Hamiltonian can 
then be made locally gauge invariant up to some power of e by introducing 
the gauge (photon) field A which transforms under the gauge transformation 
as A ^ A + Vx- It is then straight forward (see Appendix B) to show that 
the following Hamiltonian is locally gauge invariant up to C(e^): 



niQp,Qn) = H + iA.[DQ,H]--A^A^ 



= H + 6Ho{Qp,Qn) + SV{Qp,Qn), 



(4.12) 
(4.13) 



where 6Ho (SV) is the contribution from Hq (V) to the commutator terms. 
Evaluating the part SHq explicitly, we arrive at the familiar form 



'H{Qp, Qn) 



1 



N 



1=1 



(4.14) 
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where 

5V{Q,, g„) = lA ■ [Dq, V] - \nA^ \d\^, \d^q, y]] . (4.15) 

The second term on the r.h.s. of (4.12) involves the famihar form of the 
current operator in the LWL: jQ{q = 0) = -Dq]. The third term on the 
r.h.s. involves the 2-body part of the Compton scattering amplitude in the 
LWL, as will be discussed in more detail in sect. 6. 

Here we note the following point concerning the dependence of the inter- 
action term 8V on the charges and Qn- If we consider a pair of nucleons 
denoted by (1, 2), the relevant term in the dipole operator is 

Dq^2 = Q^'Vi + g(2)r2 = J2 (g + T; Ag) + ^ Ag (rf^ - rj^)) , (4.16) 

where and r are the cm. and relative coordinates of the pair (1,2), 
Tz = (r]^) + rj^)) /2, and Q = Qp+Qn, Ag = Qp-Qn- Using this expression 
in (4.15), we see that, if the potential commutes with the 2-body cm. coor- 
dinate {[V, R] — 0), the 2-body current (first term on the r.h.s. of eq.(4.15)) 
depends only on the difference of the charges Ag and not on the sum Q. 
From the Jacobi identity for double commutators it is seen that the same 
holds also for the two-body part of the photon scattering amplitude (second 
term on the r.h.s. of eq.(4.15)) ^. This is important for our later discussions, 
since the effective charges for El transitions, Qp — Cp N/A, g„ = — Z/A, 

'^This is seen as follows: First, if we apply the Jacobi identity to the double commutator 
[T^, [Dl^,V]l we obtain [T^, [D'q, V]] = 0. Second, if we apply it to [R\ [V, D^]], we obtain 
[-R% [V^, -Dq]] = ii V commutes with i?*. From these two observations, it follows that in 

the double commutator [Dq, [Dq,V]] the term proportional to the cm. coordinate in the 
dipole operator (4.16) docs not contribute, and therefore this double commutator depends 
only of the difference AQ of the charges. 
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satisfy Aq = Cp as do the physical electric charges. Therefore, if the two- 
body potential commutes with the two-body cm. coordinate, the interaction 
terms in the Hamiltonian (4.14) are invariant under the replacement of the 
physical charges are by the effective ones: SV{qp, Qn) — SV{ep, e„). 

The introduction of the effective charges is then done as usual: First, for 
the physical charges we have 

Z 2 Z 

n{ep, e„) = //o - ^ E A • ^ + ^ E + e„) . (4.17) 

Second, in order to remove the effect of the center of mass motion on the 
electromagnetic interaction, we replace the momenta in the second term 
of (4.17) by the Jacobi momenta Pj — P/A, where P = Y^iPi is the total 
momentum [37]. This replacement naturally leads to the effective charges 

9p = ^ep, Qn^-^ep. (4.18) 

Third, since A is a constant in the LWL considered here, we can use the 
following identity for the third term on the r.h.s. of (4.17): 

Zel = ^ + Zql + Nql (4.19) 

to obtain the final form of the Hamiltonian as 

niqv. Qn) = ^^' + nQp, Qn). (4.20) 
Here the first term describes the Thomson scattering of the nucleus as a 
whole [25], and 7i(gp, qn) is given by (4.14) with the effective charges (4.18). 

The procedure to remove the cm. motion therefore leads to the Hamil- 
tonian H{qp,qn), which is locally gauge invariant under the gauge transfor- 
mation involving the effective charges. Therefore one can continue to use 
all consequences of local gauge invariance, hke Ward identities etc., for the 
theory described by the Hamiltonian TC{qp, qn)- 
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4.2 The El sum rule 



The strength function (cross section) for the absorption of unpolarized pho- 
tons by a nucleus in its ground state |0) is given by 



Siu) = An'j:\{n\j{q)^\0)\'5{u;-u;r.o) (4.21) 



UJ 

LWL . O'sr^,, i-./ „^ 1 



47r'E = 0)^|0)|2 5(a; - Uno), (4.22) 



where j{q) denotes the Fourier transform of the 3-component (or any other 
space component) of the electromagnetic current operator for the effective 
charges (4.18), and ojno — En — Eq is the excitation energy of the state \n). 
The factor ^ in (4.21) and (4.22) originates from the normalization of the 
photon field. The LWL indicated in (4.22) holds if \q\R « 1, where R is 
the nuclear radius^ ^. 

In the LWL one can use current conservation to express the current op- 
erator as j{q = 0) = i[H,Dq], where Dq is the 3-component of the dipole 
operator involving the effective charges (4.18). Then one arrives at the "non- 
retarded" El sum rule [25] 

roo ^ 1 

S= dcoSiu;) = 471^5: — |(n|j(q = 0)|0)|2 (4.23) 

= -2n^0\[[H,D^],D,]\0). (4.24) 



^°We mark the current operators and quasiparticle currents for the El effective charges 
by hats in order to distinguish them from the quantities related to the physical charges. 

"'^"'^The multipole expansion of (4.21) can be performed as usual by the expansion of the 
factor exp(— • r) which defines the Fourier transform of the current operator. This leads 

to the J-multipole operators of the form [j'^^' x F^'^] j;(gr). The El multipole (J = 1) 
therefore contains, besides the LWL, also "retardation terms" oc (gi?)^, which come from 
the expansion of the Bessel functions for I = and 1 = 2. Therefore, the El multipole 
contains also terms which go beyond the LWL. 
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Of course, this sum rule is not directly observable, since it includes only the 
unretarded El multipole, and for high cu {— \q\ for a physical photon) the 
LWL is not valid. It is therefore important to identify a part of the sum rule 
which holds in the low energy region where the LWL is justified. 

We now relate the strength function and sum rule to the Fourier transform 
of the correlation function (see eq. (2.8)) 

n^-'{q',q;u;)^i T dr e^^ {0\T (f{q' ,t')f{-q,t)) |0), (4.25) 

where t — t' — t , and q and q' are the incoming and outgoing momenta. For 
the case of forward scattering we define U.^'^{q, u) = U.^'^{q, q; lo), and obtain 
from eq.(4.25) the spectral representation for the 33-component 

U{q,u;) = - El(0|i(9)l^)r (;-— - J ■ (4-26) 

„ \0J — OJnO + iO OJ + OJnO — 10/ 

Comparison with eqs. (4.21) and (4.23) gives 

An 

S{u) = — lmU{q = 0,Lj) {cu > 0) (4.27) 

(jj 

S = 2%^ U{q = 0,uj = 0). (4.28) 

So far, all relations of this subsection hold for finite nuclei. The strength 
function in the LWL, eq. (4.27), is actually meaningless in nuclear matter, 
since there is no scale which corresponds to the nuclear radius R. The sum 
rule (4.28), however, represents a bulk property of the nucleus which should 
not depend on the details of nuclear structure, and therefore the discussion 
of the sum rule in nuclear matter has been the subject of many previous 
investigations [3]- [7]. 

In order to apply cq.(4.28) in nuclear matter and to discuss the K-gi 
relation, we have to know the order in which the limits q — > 0, a; — > should 
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be taken. This is important, because if q ^ is taken first (the "a;-hmit" of 
sect. 3), the p-h excitations do not contribute to the polarization, while for 
a; — > first (the "static limit" ) they contribute at the Fermi surface [26] . We 
therefore present an alternative derivation of eq.(4.28), which makes clear 
the order of the hmits. We make use of the following identity [38]: 

^ 5^(01^ |0> = -(0| [r{x'),r{x)] |0) 5(r)) . 

(4.29) 

The use of a partial integration in r and the Heisenberg equation of motion 
for j° then gives the following Ward-Takahashi identity for the correlation 
function (4.25): 

q'^q,U^%q',q;u;) = -(0| [[i^, j°(-q)] , j°(q')] \0) + ^ (0| [?(-q),?(g')] |0) 

(4.30) 

where q''^ = {u!,q') and q^ = {uj,q)- We set a; = in this identity, then 
d d 

^PPly T77 7^ — to both sides, and finally let q' and q go to zero to obtain the 

dqi dqj 
low-energy theorem 

n(q = 0, = 0) = -(0| [[H, D,] , D^] |0), (4.31) 

where now the limit is defined as ^ first followed by g — ("static 
limit"). Comparison with eq. (4.24) then gives again the relation (4.28) with 
this particular "static limit" prescription for the low frequency / low wave 
length limit. 

Comparison with eq.(4.12) of the previous section shows that the double 
commutator in (4.31), which determines the sum rule (4.24), is the "seagull 
part" of the low energy photon scattering amplitude. This point will be 
further discussed in sect. 6. 
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4.3 The K — g£ relation 



The polarization, which determines the sum rule by (4.28), has been split 
into two parts in eq.(2.10). We will now show that the part Ha, which has 
been given explicitly in (2.21), is related to the orbital g-factors. The part 
n^, on the other hand, has no relation to the g-factors. 

Taking the static limit (a; ^ first, then g ^ 0) of Ua in (2.21), we 
obtain 

Z N 



UA{uJ = 0,q = Q)=j. 



\Jp ~ Jn 



(4.32) 



^ vf{p)pf{p) ^ " VF{n)pF{n)' 

where the currents j^'^^ and j*^^*' at the Fermi surface refer to the effective 
charges, in contrast to the currents for the physical charges used in sect. 3. 

We now follow the same steps as in sect. 3: The Ward identities determine 
the currents in the static limit as 



Relation (2.22) then gives the currents in the a;-limit: 



jp 



Jn 



Qp 



vf{p) + 



vf_ fpApY 

6 [ Pf , 



VF{n) + 



Vf_ ( PFjpY 

6 [ Pf , 



Fi{pp) 



Fi{nn) 



6 \ Pf 



(4.33) 



Vf PF[n \ TP r \ 
+ Qn— I — — 1 Fi{pn) 



(4.34) 



+ qp 



vf_ ( Pf{p) ' 
6 I _ 



Fi{pn) ■ 



(4.35) 



If we use the Landau effective mass relations of Appendix A for the nonrel- 
ativistic case, we obtain the very simple results 



Jp Jp ^ 



" AM ' 



Jn Jn ^ 



AM ' 



(4.36) 
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where the currents j^^^ and j^'^^ for the physical charges are given by eqs. 
(3.7) and (3.8) in the nonrelativistic limit fip — fJ>n — M. The relations (4.36) 
clearly show the role of the effective charges to remove the cm. motion: For 
the case N — Z the subtractions in (4.36) make sure that there are no 
contributions from isoscalar currents to the collective motion of the system. 

Inserting the static currents (4.33) and the 6i;-currents (4.36) into (4.32), 
and using the definition of the orbital g-factors (3.1), we obtain 



YlAiu = 0, = 0) = e; 



NZ 

•AM 



(geip) - ge{n)) ■ 



The sum rule (4.28) then can be expressed as 

S = (TRK) {1 + KA + Kb) = Sa + Sb, 



(4.37) 



(4.38) 



NZ 



where (TRK) — 27r e„— — is the Thomas-Reiche-Kuhn sum rule value, and 
the enhancement factor 1 + k has been split into the two parts 1 + ka and 
Kb originating from 11^ and Ub, respectively: 

1 AM 



1 + ka 



Kb 



e^NZ 
1 AM 



UA{q = 0,uj^O)^ge{p) 
ns(g = 0,cu = 0). 



[n 



(4.39) 
(4.40) 



Using (3.12), the relation (4.39) between ka and the orbital g-factors can 
also be expressed as follows: 



N 

ge{p) = ■^ + 34'^^ 

( \ ^ 



where ka is given by 



Ka 



Mvf 
3pF 



Fi{pn) 



1 - 



N - Z 
A 



(4.41) 
(4.42) 

(4.43) 
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From our discussions in sect. 3 concerning the dependence of Fi{pn) on the 
neutron excess, we see that the deviation of ha from its value for symmetric 
nuclear matter starts with {{N — Z)/A)^. 

5 Discussion 

As we pointed out in sect. 3, the expressions (3.10) and (3.11) are exact in 
nuclear matter and therefore take into account all possible contributions from 
meson exchange currents and configuration mixings, except for the conven- 
tional RPA-type contributions (e.g., the first order configuration mixing). 
The El enhancement factor /t, on the other hand, has been split into the 
pieces ka and hb in eq.(4.38), and only ka is related to the orbital g-factors. 
In this section we wish to discuss which processes are taken into account by 
kai and whether ka can be related to the experimentally measured sum rule. 

For this purpose, we return to eq.(2.10), where the polarization has been 
split into the two pieces 11^ and 11^. According to the optical theorem (4.27), 
the strength function for finite nuclei in the LWL then splits into two pieces 
as well: S{uj) = 5'a(^) + iS'b(cj). We first wish to discuss the piece 11^ in 
terms of Bethe-Goldstone diagrams. From the structure of eqs. (2.6), (2.10), 
and the energy denominator appearing in the quantity A of (2.18), it is clear 
that the piece Ua can be related only to those time-ordered diagrams which 
have a p-h cut in each internal loop integral appearing in the RPA series, 
i.e., which can be made disconnected by cutting a p-h pair in any internal 
loop of the RPA series. Diagrams without p-h cuts contribute exclusively to 
n^, see Fig. 5 for examples. 

Let us denote the contribution of a particular time-ordered diagram i, 
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which involves p-h cuts, to the total polarization by 11'^*'', see Fig. 4 for 
examples. Besides the p-h cuts, which correspond to "low energy denom- 
inators" of the form (lo — ± iS^ , the diagram i will in general also 
involve "high energy denominators" (u — E'^'sh _|_ Here E^^ denotes 

the excitation energy of an intermediate p-h pair, and E^'s^ is the excita- 
tion energy of a more comphcated intermediate state like a 2p-2h state etc. 
These high energy denominators appear in the effective vertex r^*^) (second 
row in Fig. 4) and in the effective interaction T'^^^ (third row of Fig. 4). The 
diagram i will therefore in general give contributions to both Ua and n^, 
that is, n« = + n^^. in Appendix C we show that the piece U\ can be 
obtained from the full H^*) by approximating the high energy denominators 



1. In the energy denominators which come from the effective vertex F^'^^ 
the replacement u; — > E"^^ is introduced. That is, u; is replaced by the 
excitation energy of the p-h pair which enters (or leaves) the vertex 
F^'^) in the loop under consideration. 

2. In the energy denominators which come from the effective interaction 
T^'^\ the replacement u ^{Ef^ + E^^) is introduced. That is, u; is 
replaced by the average of the excitation energies of the p-h pairs which 
enter {Ef ) and leave (Ef ) the block T^"^). 

3. If the particle or the hole in the intermediate state under consideration 
(excitation energy E"^^) has an energy dependent (dispersive) self en- 
ergy, the associated high energy denominators are replaced by the first 
two terms of their expansion around lu = E^^. That is, the effects of 
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dispersive self energies are included only in the quasiparticle energies 
and residues. 

To summarize these prescriptions: In order to obtain the contribution of a 
time-ordered graph, which has a p-h cut, to 11^, one has to replace uj in the 
high energy denominators by a value which is related to the excitation energy 
of the intermediate p-h pair. 

Let us then consider the contribution of a particular time-ordered graph, 
n*^'), which has both p-h and higher energy cuts, to the LWL sum rule. In 
Appendix C we show from analyticity that the following relation^^ holds for 
any diagram i: 

5« = 47r /°° — ImnW(q = 0,a;) = 27r2n«(q = 0,u; = 0). (5.1) 

Since the analyticity arguments leading to (5.1) are not invalidated by the 
replacement of the high energy denominators by a;-independent quantities 
following the hues discussed above, a relation similar to (5.1) holds for 11^^ 
(and therefore also for II^'') separately: 

=47r /°°— Imn5^(q = 0,u;) = 27r2n?(q = 0,cu = 0). (5.2) 
Jo UJ 

On the other hand, by using the above prescriptions to obtain H^-* from the 
full n('\ we see that, as long as the energies E'^'s^ are on the average large 
compared to the typical p-h excitation energies E'p^, the following relation 
holds: 

n§^(g = o,cj = 0) ~ n(^)(q = o,cj = 0). (5.3) 

^^As long as we consider a particular diagram z, the imaginary part Imll'^*) is not 
necessarily positive. 
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Here the symbol ~ in (5.3) indicates that this relation is approximately valid 
in finite nuclei as long as in the average 

£;high ^ph^ ^Yiile in nuclear 
matter it is exact because for q = we have = from momentum 
conservation. Combining with eq.(5.1) we obtain the following relation for 
the contribution of the diagram i to the sum rule: 

5« ~ Sf. (5.4) 

This relation indicates that the higher excited states (2p-2h etc.), which are 
mixed into the low energy p-h states in the diagram i under consideration, 
contribute to the sum rule mainly via their real parts. Their imaginary 
parts lead to deviations of the sum rule 5"^*^ from the value S^2\ but these 
deviations are small, although the strength function (S^''\uj)) itself is, off 
course, influenced by these imaginary parts, which give rise to the well known 
"spreading widths" . 

Summing over all diagrams i which have p-h cuts and possibly also higher 
energy cuts, we then have the following relation for their contribution to the 
sum rule: 

^ 5» c^Sa^ (TRK) {ge{p) - ge{n)) , (5.5) 

i 

Our above discussion indicates that, although in principle the k — ge relation 
takes into account the mixing of the higher excited states (2p-2h etc.) into the 
p-h states only via their real parts, the contribution of the imaginary parts of 
the higher excited states to the sum rule is comparatively small. Therefore 
we can say that, besides the contributions from the p-h cuts, the k — g£ 
relation includes also the most important effects of the mixing between the 
p-h and the higher excited states. In the same approximation, the remaining 
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piece Kb arises exclusively from those diagrams which have no p-h cuts at 
all, like those shown in Fig. 5. 

Since the part ka includes the p-h cuts and the most important part of the 
coupling to higher excited states, one can expect that ka will be dominated by 
the GDR, which is a superposition of collective p-h pairs. We can understand 
this from eq.(2.10), which shows that the part Ua involves the total vertex 
r, which is obtained from the effective p-h vertex T^^^ via the RPA-type 
equation (2.6). The vertex F is thus generated by a collective superposition 
of p-h states represented by the p-h propagator A, and contains the discrete 
poles as well as the continuum cut due to the collective p-h pairs. The discrete 
poles correspond to the "zero sound" modes in nuclear matter, and to the 
giant resonances in finite nuclei. Therefore, the piece ka includes the effect 
of the most prominent excitations, that is the collective p-h excitations, in 
the low energy region, and also the most important part of their mixing to 
the higher excited states. 

On the other hand, the piece kb arises from those diagrams which have 
no p-h cuts, like the examples shown in Fig. 5. It is known that, because of 
the short range nature of the tensor force between the nucleons, the "effective 
energy denominators" for these diagrams are large, typically several hundred 
MeV [6, 30]. These diagrams will therefore contribute to the strength func- 
tion mainly in the region well beyond the GDR. In this high energy region, 
the LWL is no longer valid, that is, the contribution of these diagrams to the 
LWL sum rule (4.28) has no connection to the measured strength function. 
The K — ge relation (4.39) therefore connects the "observable" part of the en- 
hancement factor in the LWL sum rule, which is 1 -|- ka, to other observable 
quantities, namely the angular momentum g-f actors. 
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The important point of our above discussions was the observation that the 
part includes the effects which come from the collective p-h excitations. 
In Appendix D we illustrate this point more explicitly for the case of nuclear 
matter. There we show that the RPA-type equation for the vertex (2.6) is 
equivalent to the Landau equation [13, 16, 28] in the vector-isovector channel, 
and derive the expression for 11^ in terms of the solutions to the Landau 
equation, see eq.(D.13). The discussions in Appendix D clearly demonstrate 
that the effects of the collective excitations of the system are included in the 
part n^. 

Let us come back to the k — ge relation (4.39), and discuss the nuclei in 
the ^°^Pb region for the sake of illustration. The theoretical values of ge{p) 
and ge{n) for nuclei in the lead region given in table 7.12 of ref. [27] are 
gi{p) ~ 1 + 0.13, ge{n) ~ —0.07. These values are consistent with (3.12), 
and with the empirical ones of ref. [4]. From Eq.(4.39) one then obtains the 
estimate ka — 0.2. In the analysis of ref. [11], which uses the experimentally 
measured scattering cross section to extract the total photoabsorption cross 
section via dispersion relations, it was concluded that "any reasonable pre- 
scription gives (experimental) values of kgdr between 0.2 and 0.3", where 
'^GDR was extracted from the area under a Lorentzian curve fitted to the 
GDR. This would indicate at least a quahtative consistency between theory 
and experiment, since our k,a can be identified with Kqdr as discussed above. 
However, the analysis of Ref. [8] indicates that kqbr might be enhanced by 
retardation corrections. This would indicate a discrepancy between theory 
and experiment, which should be further investigated 

^^Since in the discussion of refs.[8] and [11] (and also in other papers) it was assumed the 
quantity ge, which enters in the k — gg relation, includes only the meson exchange current 
corrections to the orbital g-factor and not the other nuclear structure effects (configuration 
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6 Connection to the photon scattering am- 
plitude 



In this section we wish to discuss the physical reason which underhes the 
K — g£ relation. In the original discussions on this relation [3]-[7], the en- 
hancement factor K was calculated in perturbation theory directly from the 
commutator of eq.(4.24), and the orbital g-factor was calculated from the 
current according to eq.(3.1) in the same order of perturbation theory. In 
this approach, the k — relation appears a bit fortuitously. On the other 
hand, following the lines discussed in sect. 4, it is clear that the low energy 
theorem (4.31) directly leads to the k ^ gc relation: The polarization, which 
appears on the l.h.s. of this relation, contains the part 11^, which is related to 
the p-h excitations as illustrated in Fig. 3, and in the hmit of a; — > followed 
by \q\ — ^ these p-h excitations can take place only at the Fermi surface. 
It is then immediately clear from Fig. 3 that this process is proportional to 
the matrix element of the effective vertex F'^'^^ (the black square in Fig. 3) 
at the Fermi surface, which directly gives the quasiparticle current and the 
orbital g-factors via the definitions (2.19) and (3.1). On the other hand, the 
r.h.s. of (4.31) is related to k, which establishes the desired relation. 

From our discussions in subsect. 4.2 it is clear that the basic relation 
(4.31) follows from gauge invariance. Actually it can be viewed as a con- 
mixings etc), it was concluded that theory would be consistent with the larger value for 
KGDR- (For example, the meson exchange current corrections in table 7.12 of ref.[27] 
are on the average Sgi{p)'^^^°'^ ~ 0.28 and (n)™®^°" ~ —0.16, and also the results of 
ref. [29] are very similar. If eq.(4.39) would hold only for the mesonic corrections to ge, 
one would obtain a large value ka — 0.44.) However, as is clear from our derivation, 
the orbital g-factors which enter in the k — gi relation are the total ones, including the 
configuration mixing effects [36] (except for the conventional RPA-type effects like first 
order configuration mixing), besides meson exchange currents. 
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sequence of current conservation applied to the nuclear Compton scattering 
amplitude: Eq.(4.20) shows that the Compton amplitude naturally splits into 
the Thomson amphtude, which follows from the first term in (4.20), and the 
"intrinsic" scattering amplitude, which follows from the Hamiltonian (4.14) 
by using the effective charges (4.18). This "intrinsic" scattering amphtude 
7]nt consists of two pieces. 



where TL^'^ is the current-current correlation function defined in (2.8) and is 
often called the "resonance amplitude" in this connection [8], and Mj^ is 
the "seagull" amplitude described by the double commutator term in (4.12). 
Some examples for are shown in Fig. 6. 

Since we have shown that the underlying Hamiltonian (4.14) is invari- 
ant with respect to local gauge transformations generated by the effective 
charges, the amplitude (6.6) must satisfy current conservation, in particular 



Since u and q in this relation can be treated as independent quantities, the 
limit cu ^ first followed by q ^ leads to 



Ti^iq = 0,00 = 0)= W^{q = 0,u; = 0) + Mg{q = 0,u; = 0) = . 

(6.8) 

Since M'^i{q = 0,lu = 0) = {0\[[H, D% D^O) follows from (4.12), we imme- 
diately obtain low energy theorem (4.31). The underlying physical reason 



^''The current operators which define the correlation function n'^" include the 2-body 
exchange currents as specified by the single commutator term in 5V of eq. (4.15). The 
piece Mf^ therefore contains all processes which cannot be expressed as the product of 
two current operators. 



(6.6) 



q^.quT-Z{q.uj)^0. 



(6.7) 
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Figure 6: Some selected examples for diagrams contributing to the seagull 
amplitude Mint of eq.(6.6). The seagull-type vertices in the first and third 
diagrams come from the terms oc in the non- interacting part of the Hamil- 
tonian (4.14), while the vertices in the second and fourth diagrams come from 
the interaction part 8V . 

for the K — gi relation can therefore be identified as the condition of current 
conservation for the nuclear Compton scattering amplitude. 

In Appendix E we discuss further connections between the low energy 
Compton scattering amplitude and the orbital g-factors, making contact to 
the analysis in ref.[8] of photon scattering data. There we also comment on 
the sum rule derived a long time ago by Gerasimov [31]. It is well known 
that this relation seemed to indicate that the (unrctardcd) El sum rule has 
the same value as the integral over the total photoabsorption cross section 
including retardation effects and all higher order multipoles (see Appendix 
E). However, its vahdity has been questioned in a series of papers [33]-[35]. 
In particular, in ref. [35] it has been shown that the assumption concerning 
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the high energy behavior of the amphtude is actually not valid due to the 
presence of the anomalous magnetic moment term. 

7 Summary and conclusions 

In this work we used the Landau-Migdal theory to discuss the orbital angular 
momentum g-factor of a quasiparticle and the El sum rule for isospin asym- 
metric nuclear matter. The relations obtained for the orbital g-factors are 
in principle exact and hold also in relativistic field theory. For the El sum 
rule, we had to restrict ourselves to a nonrelativistic framework because of 
the problems arising from the center of mass motion. We discussed in detail 
the form of the Hamiltonian which is invariant with respect to local gauge 
transformations generated by the effective El charges. 

We have split the strength function into two parts, where one comes from 
the p-h cuts including the effects of the higher excited states via their real 
parts, and the other comes from cuts at higher excitation energies. We have 
shown generally that the former part is related to the orbital g-factors, while 
the latter part has no relation to them. The former part has a close relation to 
the collective excitations of the system, i.e., the zero sound modes in infinite 
systems and the giant resonances in finite nuclei. We have discussed the 
importance of the k — relation, which effectively separates the observable 
part of the LWL sum rule, which is related to the strength function in the 
low energy region, from the rest. Our discussions, which do not rely on 
perturbation theory, can serve to put many previous investigations on the 
K — ge relation on a theoretically firm basis. 

Concerning possible extensions, we would like to remark the following 
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points: First, the methods used here to relate k to ge refer to infinite nuclear 
matter, and it would be interesting to investigate to what extent they can be 
applied also to finite nuclei. Second, as we mentioned in the Introduction, 
very interesting attempts are now being made to extend the range of applica- 
bility of the Landau-Migdal theory [24] to give a more general description of 
nuclear collective vibrations. The basic idea is to generafize the definition of 
the quantity A, which appears in the equation of the vertex (2.6) etc., so as 
to include also more complicated configurations. It would be very interesting 
to see whether the results derived in this paper can be extended according 
to these lines. 
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Appendices 



A Ward identities and the relativist ic Lan- 
dau relations for N ^ Z 

In this Appendix we derive some relations used in sect. 3, which are based on 
gauge and Lorentz invariance. For more detailed discussions for the case of 
isospin symmetric nuclear matter, we refer to refs.[26, 40]. 

A.l Ward identities and currents 

If Qp and Qn are the bare electric charges for protons and neutrons, the 
requirement of local gauge invariance leads to the Ward-Takahashi identity 
between the electromagnetic vertex Fq, and the propagator Sa for protons 
{a — p) and neutrons {a — n): 

q,K{k\ k) = {s-\k') - S-\k)) . (A.l) 

By setting — {) first and then letting q — > 0, one obtains the Ward identity 
for the vertex in the "static limit" [13] : 

Tt'\Kk)^-Q^VkS-\k). (A.2) 

If the propagator has a quasiparticle pole ko — ea{k) with residue Za{k), we 

can expand S^^{k) = Za{k)^^ {ko — ea{k)) + O {ko ~ ea{k))^ to obtain the 
current in the static limit on the quasiparticle energy shell (/cq — ^ ^a{k)) as 
foUows^^: 

Jf\k) = Z^{k) T^:'\k, k) = Vk ea{k) (A.3) 



^^In this Appendix, we will denote the current for the charges Qa as J, which cor- 
responds to j of sect. 3 for the physical charges, and to j of sect. 4 for the effective 
charges. 
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On the Fermi surface, this gives the static current as J^^^\k) = hQaVpicx), 
where VF{a) is the Fermi velocity for protons or neutrons. 

If we take the static hmit {qo — > first, then |qf| — > 0) in eq.(2.22) for 
|fc| on the Fermi surface, and use the form (3.3) for the p-h propagator, we 
can use the above form of the current in the static hmit on the Lh.s. and 
the r.h.s. under the integral. In this kinematics the interaction t^^^ depends 
only on the angle between the directions of k and £. One then obtains the 
current in the "a; limit" as 




where / is the spin independent part of t^^\ The angular integrations are 
trivial, and if one uses the definition of the dimensionless interaction (3.4) 
one obtains the currents (3.5) and (3.6) for the physical charges, or (4.34) 
and (4.35) for the effective charges, where Fi are the coefficient of the i — 1 
Legendre polynomial in the expansion of the spin independent part of the 
interaction 

A. 2 Lorentz (Galilei) invariance for N ^ Z 

The requirements of Lorentz invariance have been discussed in ref.[40], and 
are easily generalized to the case N Z as follows: 

If we observe matter from a system S which moves with velocity —v rela- 
tive to the system S where matter is at rest, the Fermi distribution function 
of neutrons (a = n) or protons {a = p) for fixed momentum £ will change 
by an amount n^(«) — Ti£{a), where n and n are the Fermi distributions in 
S and S. According to Landau's basic hypothesis [15], this corresponds to a 
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change of the quasiparticle energy (for fixed momentum k) by an amount^^ 



ik'ia) - ek'{a) = 2 



(27r)^ 



Y: fap{k\£){he{P)-ne{P)) , (A.l) 



I3=p,n 



where fa/s is the spin independent part of the quasiparticle interaction 

as in (A. 4). On the other hand, the l.h.s. of (A.l) must agree with the result 

obtained from a Lorentz transformation: 

ik'ia) ^ -f {ek{a) + k ■ v) , (A.2) 

where ^ — {1 — v"^)^ , and k' and k are related by the Lorentz transformation 

fc' = fe + ek{a)v-f - v (t) ■ fe) (1 - 7) . (A. 3) 

For the r.h.s. of (A.l), we use the fact that the Fermi distribution is Lorentz 
invariant: 



ne'la = nela 



(A.4) 



where £' and £ are related by a Lorentz transformation analogous to (A.3). 
To first order in v one easily obtains from (A.2), (A.3) and (A.4): 



ik'ia) - ek'{a) = k ■ v - ek{a)v ■ Vkia) 

ni{a)-ni{a) = ei{a)v ■ £ 6{\i\ - pria)) , 



(A.5) 
(A.6) 



where Vkicx) — ^k^k{(^) is the quasiparticle velocity. Inserting the relations 
(A.5) and (A.6) into (A.l) and setting |fc| —pF{a), we obtain 



Vf{p) + y 



. Pf . 



Flipp) + ^ 



fJ'n ( PF{n)' 



fip \ Pf 



Fiipn) 



Pf{p) 



^^Since the quasiparticle interaction fai3{k',i) used in (A.l) refers to the system S, 
the relation (A.l) actually holds only up to the first order in v, which is sufficient for 
our purpose. Concerning the Lorentz transformation of the quasiparticle interaction, see 
ref.[40]. 
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VF{n) + — 
b 



K Pf J I^uKpfJ 



PF[n) 

IJ'n 

(A.7) 



where /Xq, = eF{a) are the chemical potentials (Fermi energies) of protons 
and neutrons. These relations were used to derive eqs.(3.7), (3.8) and (4.36) 
in the main text. 

B Local gauge invariance of the Hamiltonian 
(4.12) 

Here^^ we show the invariance of the Hamiltonian (4.12) under local gauge 
transformations up to C(e^). For this, we note that the arguments given 
below eq.(4.8) in the main text show that for any operator B, which commutes 
with the z-component of the total isospin operator {[Tz,B] = 0), the gauge 
transformations like (4.7) or (4.8) can be expressed by the dipole operator. 
That is, if [T^, B] = 0, the gauge transformation of B can be expressed as 
follows: 

B exp [-ic ■ D] Bexp [ic-D]=B- id[D\ B] - -dd[D\ [D\ B]] + ... 

(B.l) 

We have to apply this transformation, together with 

A^A + c (B.2) 

to the Hamiltonian (4.12). First, since \Tz.,H] = 0, we can use eq.(B.l) for 
B — H, as we have done already in (4.11) of the main text. Second, also for 

-^''To simplify the notations of this Appendix, we will not indicate the dependence of the 
Hamiltonian and the dipole operator on the charges Qp and Q„. 
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the single commutator term on the r.h.s. of (4.12) we have [T^, [D^, H]] = 0, 
as can be seen easily from the Jacobi identity for double commutators. Third, 
the last term in (4.12) is already of 0{e^), and therefore it is affected only by 
the gauge transformation of A, eq.(B.2). We therefore obtain for the gauge 
transformation of the Hamiltonian (4.12): 

n exp [-ic ■D]{H + i (A^ + c^) [D\ exp [ic ■ D] 

- ^ (^M^' + 2AV + cV) [D\ [D\ H]] (B.3) 

Using then (B.l) for B — H and B — i[D\H], it is easy to see that we 
simply get back the three terms on the r.h.s. of (4.12), which shows the 
gauge invariance of 7i. 

C Time-ordered diagrams 

In this Appendix, we first explain the three points stated in sect. 5 concern- 
ing the relation between H^, which was originally defined in the language 
of Feynman diagrams in sect. 2, and the time-ordered Bethe-Goldstone dia- 
grams. Although the arguments presented below can be generalized by using 
dispersion representations for the vertex, the interaction and the self energy, 
we prefer to discuss representative examples in order to be definite. The ar- 
guments can be applied immediately also to other cases. We then will derive 
the relation (5.2). 

C.l 2p-2h states induced by vertex corrections 

The third diagram in Fig. 4 arises from the contribution to the quasiparticle 
vertex F^'^^ which is shown in Fig. 7. By performing the frequency integrals 
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for the two loops, one obtains the energy denominators 
where the energies Ei and E2 are indicated in Fig. 7. 



(C.l) 




Ej E2+C0 



Figure 7: A contribution to the vertex r^^\ which is associated with the 
factor (C.l). 



If this vertex appears in a Fcynman diagram together with the p-h prop- 
agator A of eq.(2.18), the effect of A is to replace the energy denominator 
(C.l) by 



(e;,^ - E,) (ex_ - E2) 



(C.2) 



This is the "high energy denominator" of IIa. On the other hand, the high 
energy denominator in the time-ordered diagram (third diagram of Fig. 4) 
is 



(uj + ex_ - El) (ex_ - ^2) 



(C.3) 
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If we replace uj E^^ = e\_^ — e\_ in (C.3), the expression agrees with 
(C.2). In other words, if we replace in the high energy denominator of the 
time-ordered graph cu by the excitation energy of the p-h pair which enters 
(or leaves) the vertex, we obtain the contribution of this graph to 11^. This 
is the content of point 1. discussed in sect. 5. 



C.2 2p-2h states induced by the interaction 

The fifth diagram in Fig. 4 arises from the contribution to the quasiparticle 
interaction T^^^ which is shown in Fig. 8. By performing the frequency 
integral for the loop, one obtains the energy denominator 



[ko-io-E]-' , 



(C.4) 



where the energy E is indicated in Fig. 8. 




( » ko+ ) 



( ko - 



Figure 8: A contribution to the interaction T'^^\ which is associated with the 
factor (C.4). 



If this interaction appears in a Feynman diagram together with two p-h 
propagators A, which correspond to the incoming and outgoing p-h pairs. 
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the effect of ^4 (see eq.(2.18)) is to replace tlie energy denominator (C.l) by 



^ (eA+ +eA_ -ca; -eA'_) -E 



-1 



2, -J - (C-8) 

Tliis is tlie "liigli energy denominator" of H^. On the other hand, the high 
energy denominator in the time-ordered diagram (fifth diagram of Fig. 4) is 

'uj + ex_-ey^-E[^ . (C.6) 

If we replace u ^ \ (sf + Ef) = \ (eA+ - eA_ + ca; - eA'_) in (C.6), the 
expression agrees with (C.5). In other words, if we replace in the high energy 
denominator of the time-ordered graph uj by the average of the excitation 
energies of the p-h pair which enters and leaves the interaction, wc obtain 
the contribution of this graph to 11^. This is the content of point 2. discussed 
in sect. 5. 

C.3 2p-2h states induced by the self energy 

The p-h propagator A defined in sect. 2 includes by definition the effects of the 

self energies on the quasiparticlc energies and the residues. For example, if 
wc denote by 7 the non-interacting electromagnetic vertex, the contribution 
(^7/17) to Ha in eq.(2.10) includes the following factor which is associated 
with the diagram shown in Fig. 9a: 

s(£;,j ^ ^ as/afeo(£;Aj _ 



(u;-£;a++£;a_) ^-^a++£^a 

Here the energies E\ do not include the effect of the dispersive self energy, 
which is shown in Fig. 9b. After performing the frequency integrals for the 
two loops in the Feynman diagram for S, one obtains the energy denominator 



ko + --E) , (C.8) 
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where the energy E is indicated in Fig. 9b. 





E 



(a) 



(b) 



Figure 9: (a) A contribution of dispersive self energies to the polarization, 
(b) A dispersive self energy part, which is associated with the factor (C.8). 

Using this in eq.(C.7), we see that the high energy denominators in Ua 
have the form 



On the other hand, the high energy denominators of the time-order diagram 
for the polarization are 



Expanding the high energy denominator (second factor in (C.IO)) around uj = 
— E\_ and taking only the first two terms of this expansion, we obtain 
an expression which agrees with (C.9). In other words, if we replace the high 
energy denominator of the time-ordered graph by its first two terms of an 
expansion around the excitation energy of the p-h pair under consideration. 




(C.9) 




(C.IO) 
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we obtain the contribution of this graph to Ha- This is the content of point 
3. discussed in sect. 5. 



C.4 Proof of relation (5.1) 

For each time-order graph which contributes to the polarization, there is 
another one which is obtained from the original graph by crossing the pho- 
ton lines. That is, each graph n^*^ which is characterized by a particular 
topology of the internal propagators and two-body interactions, can be split 
into two parts, n^'^ = n+'' -|- n^-*, where U^^ does not involve crossed pho- 
tons ( "forward term" ) , and 11^^ involves crossed photons ( "backward term" ) . 
By definition, the intermediate states in H^^ have either no photons (energy 
denominators of the form {cu — En + ie)~^) , or one photon (energy denomina- 
tors are independent of u) . On the other hand, the intermediate states in H^' 
have cither two photons (energy denominators of the form {u + En — ie)~^), 
or one photon (energy denominators are independent of a;). It follows that 
n^^(a;) is analytic in the upper oo plane [39], and has a non- vanishing imag- 
inary part only for real positive cu. Prom these properties it follows that 
n+^(a;) satisfies a dispersion relation 

n»(a;) = ^ f . Im n«(a;0 . (C.ll) 

Since n^''(a;) and n^'*(a;) differ only by the external photon lines, they give 
the same contribution at a; = 0, that is, = 0) = n^''(a; = 0). We 

therefore obtain from eq.(C.ll) in the limit u — 

r^Imn«(a;)= r^Imn«(c.) 

Jo CU Jo CO 

= 7rnJ^(u; = 0) = ^n«(a; = 0), (C.12) 
50 



which is the content of eq.(5.1). 



D Connection between Ua and the collective 
excitations 

In this Appendix we derive the expression for the polarization in terms of 
the solutions to the Landau equation, referring for simplicity to symmetric 
nuclear matter (N=Z) and low energy and momentum transfer, where 11^ of 
eq.(2.21) depends only on the ratio s — Lij/{\q\vF)- For the case N=Z, the 
effective charges (4.18) become ±|ep for protons and neutrons, and therefore 
the currents in the polarization (2.21) are the isovector (IV) currents. 
If we approximate the p-h propagator in (2.21) by (3.3), the quantity 

OJ - q ■ Vk + 

(D.l) 

satisfies the equation 
{s - cos Ofc) (/>(fc, s) = jPik) + / dQe F'{k ■ 1) 0(^, s) , 

(D.2) 

where cos Qk = Q ■ k, F'{k ■ £) is the spin independent part oc (ti • T2) of the 
dimensionless Landau-Migdal interaction (3.4), and j'iv (^) isovector 
part of the quasiparticle current in the a;-limit at the Fermi surface, i.e., from 
eqs. (3.7), (3.8): 

ji-\k) ^ kji^^ ^kl[^-'fiF,- Fi)) . (D.3) 

For small \q\ it is sufficient to seek the solution to (D.2) in the form 0(fc, s) = 
k(j){k,s), where the scalar function (j){k,s) satisfies an equation similar to 
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(D.2), but with the interaction retarded by one unit of angular momentum: 
T{k ■ £) —>■ T{k ■ £) k ■ i. We can write the equation for s) in the form 

/ dQ, K{k, I) </^{£, s) + S^^j-) = s (j>{k s) (D.4) 

with the kernel 

K{k, I) = cos Gfe 5{nk - n^) + ^^^^^F'{k ■e)k-l (D.5) 

One can express the solution to the inhomogeneous equation (D.4) by the 
solutions to the homogenoues equations 

J dneK{k,£)^£,s) = s$(fe,s) (D.6) 

J dnk^{k,s)K{k,£) = s'^{£,s), (D.7) 

where we must distinguish between the right and left eigenfunctions because 

of the nonhermiticity of the kernel K. The left eigenf unction is related to 

the ordinary h.c. of $ by $(fe, s) = $'(fc,s) where Ng is chosen to 

cosGfe 

satisfy the normalization (D.8). 

For \s\ < 1 there exist "free" solutions to eq.(D.6) and (D.7), namely 
$(fe, s) = $(A;, s) = . — S{s — cos Q^), and in this case the homogeneous 
equations can be solved iteratively for any value |s| < 1, i.e., there is a con- 
tinuum of solutions. For |s| > 1, on the other hand, no free solutions exist, 
and the eigenvalue s becomes discrete. These are the zero sound solutions, 
which have been investigated in detail in many works [28, 41, 16]. 

We can choose the right and left eigenfunctions so as to form a complete 
orthonormal set: 

ldQk^k:s')Hk,s)^Ss',s (D.8) 
s) m s) = 6 {ilk- He) , (D.9) 
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where 6s'.s stands for a delta function in the case of continous eigenvalues, 
and the Kronecker delta symbol in the case of discrete ones. The symbol J2g 
refers to an integral over continuous eigenvalues, and a sum over discrete ones. 
The solution to the inhomogeneous equation (D.4) can then be expressed as 
follows: 

E.,/dn /f:';!^y' cose,. (d.io) 

where cos Qe — q- l. Returning to the correlation function (2.21), we obtain 
in the hmit of low (a;, q) the form n^(c<;, q) — S'^^IIa{uj, q), where 11^ is given 

by 

(D.U) 

where 

Ag, = |dQfe$(fc,s') 

Eg, = J dnk'^{k,s') cosQk = A*,Ns, . (D.12) 

Prom the homogeneous equations (D.6) and (D.7) it is easy to see that by 
reversing the direction of the momentum one obtains a solution with the op- 
posite sign of the eigenvalue: $(— fc, s) = $(fc, —s) and $(— fc, s) = $(fc, — s). 
One can therefore restrict the summation over the eigenvalues in (D.ll) to 
positive ones and write 

nA(^,g) = -j^ (jivT E., MN, ( ^— ■ 

AppVp ^ ' * \s — s'-|-z7 s-\-s'—i'y) 

(D.13) 
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E Some relations for the Compton scattering 
amplitude 



In this Appendix we discuss relations between the low energy Compton scat- 
tering amplitude and the orbital g- factors (or the enhancement factor 1+ka), 
and make contact to the sum rule derived a long time ago by Gerasimov [31]. 

E.l Relations for the low energy scattering amplitudes 

If we add the Thomson amplitude 

to the intrinsic part (6.6), we get the total photon scattering amplitude 

T'{q,uj) = W^{q,uj)+M^^{q,uj) (E.2) 

with M = Mint + Mtyi- Similar to the correlation function 11, we can split 
also the seagull amplitude into two parts as M = Ma + Mb- (Examples for 
Ma are shown by the first three diagrams in Fig. 6, while the fourth diagram 
contributes to Mg.) 

According to eqs.(4.37), (3.12), and (6.8), the "A-parts" of the scattering 
amplitude (E.2) can be expressed in terms of the orbital g-factors, or in terms 
of follows: 

N 7 N 7 

n^(g = 0,u; = 0) = 5^Ul—-^g,^p)-g,^n))^el^{l + KA) 

(E.3) 

M^^(g = 0,w = 0) = -n^(q = 0,a; = 0) - 



p 



MA 

f«(.)^-^(:.^). (E.4) 



54 



These relations between the low energy nuclear Compton scattering ampli- 
tudes and the El enhancement factor have been used extensively in the anal- 
ysis of experimental data [8]. One should keep in mind, however, that they 
are valid only for the "A-parts" . The "B-parts" of 11 and M, which cancel 
each other in the low energy limit, have no relations to the orbital g-factors. 

E.2 Relation to Gerasimov's sum rule 

In sect. 4 we made use of a spectral representation of the correlation function 
for fixed q, see eq.(4.26). Although for a real photon one has \q\ — OJ, in the 
energy region where uj R «1 one can approximate |qf| ~ 0, leaving uj finite. 
This leads to the unretarded El sum rule (4.23), (4.24). On the other hand, 
there exists a subtracted dispersion relation for the full forward scattering 
amplitude T{ijj) = Tij{q,u!) CiCj for a physical photon {\q\ = uj) [32]: 



This dispersion relation is often used to extract the total photoabsorption 
cross section Ctot from the measured elastic cross section [10, 11, 8]. 

According to eq.(E.2), the amplitude T consists of the "resonance ampli- 
tude" n(a;) = njj(a;) CiCj and the "seagull amphtude" M{u;) = Mij{uj) titj. 
In the energy region beyond the GDR and below the pion production thresh- 
old (o^GDR << ^ « <^7r); the energy dependence of the seagull amplitude 
M (fig. 6) is weak: M{uj) ~ M{uj = 0). Therefore one can replace T — > 11 on 
the l.h.s. of (E.5) to get the relation 



which is valid for a;GDR « io « io^n- If one makes the further assumption 
that in this energy region (a; >> ujgdr) the resonance amplitude n(a;) is 




(E.5) 




(E.6) 
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already small because of large energy denominators, one might use (E.6) 
formally also for a; — > oo to obtain 



This conjecture due to Gerasimov [31] has often been interpreted as the "can- 
cellation between retardation effects and higher order multipoles" [3],[33]- 
[35], since the l.h.s. of (E.8) involves only the unretarded El multipole, while 
the r.h.s. involves the total cross section. It has, however, been shown explic- 
itly in ref. [35] that the assumption (E.7) on the high energy behaviour of the 
resonance amplitude is not valid due to the contribution of the anomalous 
magnetic moment term to the dispersion integral. Therefore it was pointed 
out that eqs. (E.7) and (E.8) are actually not vahd. 




(E.7) 



and comparison with the exact relation (4.28) would lead to 




(E.8) 
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